Fast generation of multiparticle entangled state for flux qubits in a circle array of transmission line 

resonators with tunable coupling 
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We study a one-step approach to the fast generation of Greenberger-Home-Zeihnger (GHZ) states in a circuit 
QED system with superconducting flux qubits. The GHZ state can be generated in about 10 ns, which is much 
shorter than the coherence time of flux qubits and comparable with the time of single-qubit operation. In our 
proposal, a time-dependent microwave field is applied to a superconducting transmission line resonator (TLR) 
and displaces the resonator in a controlled manner, thus inducing indirect qubit-qubit coupling without residual 
entanglement between the qubits and the resonator. The design of a tunably coupled TLR circle array provides 
us with the potential for extending this one-step scheme to the case of many qubits coupled via several TLRs. 

PACS numbers: 03.67.Lx, 42.50.Dv, 03.67.Bg, 85.25.Cp 
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I. INTRODUCTION 

Entanglement lies at the heart of quantum mechanics and 
plays a key role in quantum information processing. In 
quantum error correction, quantum teleportation and quantum 
cryptography, the generation of high fidelity multi-particle 
entangled states, such as the Greenberger-Horne-Zeilinger 
(GHZ) state or two qubit Bell state, is required'. Therefore, 
the preparation and verification of the entangled states are of 
great practical importance in quantum information processing 
systems. 

Superconducting qubits2i^ are the most promising candi- 
date for realizing solid-state quantum information process- 
ing. Generation of the GHZ state in superconducting quantum 
circuits is consequently a highly important issue. While 14- 
particle and 10-particle entanglement have been experimen- 
tally demonstrated in trapped-ion systems^, and photonic sys- 
tems^, respectively, theoretical studies have focused on gen- 
erating the three-qubit GHZ state in superconducting charge*, 
flux^ and phase** qubit circuits with direct qubit-qubit interac- 
tion. Experimental demonstrations of entanglement have also 
been limited to the cases of the twc^"'^ or three'-'"'^ particles 
so far. Therefore, how to generate multiqubit GHZ states in 
superconducting quantum circuits is still an open question. 

The circuit QED system'^ provides a possibly scalable 
method of realizing quantum information processing with su- 
perconducting qubits. In such a system, the quantized mi- 
crowave field can act as a data bus to transfer information be- 
tween qubits. The generation of a multiparticle (three or more 
particles) GHZ state using the system of the circuit QED has 
been proposed in Refs.l l9h 24t However, the generation of the 
GHZ state in Refs. I20l - t22l is based on measurement, and is 
thus probabilistic. The probability of generation of the GHZ 
state exponentially decreases with the number of qubits. The 
proposal in Ref. [24, which is similar to that for trapped-ion 
systems^^'^* and atomic systems^^-^^, comprises deterministic 
generation of the GHZ state in a system of either flux qubits 
or charge qubits. However, it is difficult to realize this ap- 
proach in current experiments due to the existence of a few 



practical problems, as follows. First, to prepare a high-fidelity 
GHZ state, the time control of the dc current pulse should be 
precisely set around 2tx jujr, where cj,. is the frequency of the 
fundamental cavity mode and the usual range is from 27r x 1 
GHz to 27r X 10 GHz. It is a major technical challenge to carry 
out the proposed experiment with the present time resolution 
of high-performance commercial arbitrary waveform genera- 
tor, which is generally about 1 ns. Furthermore, changes in the 
half of the flux quantum in an a loop (defined in Sec.|II| with 
a typical area, for example, 10 ^in? also represent a technical 
problem, because it is necessary to apply a biased current up 
to 1 mA via on-chip biased line in an ultra-low temperature 
environment if the mutual inductance between the a loop and 
the biased line is ^ 1 pH. Second, the preparation time is one 
or even two orders of magnitude longer than that of the single- 
qubit operation. Faster preparation requires reduction of the 
frequency of the fundamental cavity mode. This may lead 
to more operational errors arising from the thermal excita- 
tions in the cavity. Third, the number of superconducting flux 
qubits that can be placed around current anti-node point of the 
one-dimensional superconducting transmission line resonator 
(TLR) is limited, and thus the proposed scalable method for 
many qubits remains an open question. 

To overcome the problems encountered in previous studies 
(e.g., in Ref. 24) and make the proposal experimentally more 
feasible, we introduce a one-step multi-qubit GHZ state gen- 
eration method in the system of a circuit QED with flux qubits. 
The advantages of our proposal are as follows. 1) The clas- 
sical driving field is directly applied to the TLR, in contrast 
to the usual case in which the driving field is separately ap- 
plied to the qubits. Thus, the interactions between the qubits 
and the classical field are induced by coherently displacing 
the cavity field. This facilitates experiments not only to obtain 
homogenous coupling constants, but also to realize synchro- 
nization of the driving on different qubits. 2) The generation 
time can be as short as the single-qubit operation time. 3) Our 
proposal can be extended a case of many qubits case by using 
a coupled TLR array as a data bus. We expect that our pro- 
posal will work well for more than 20 qubits from our current 
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experiments with the simplest setup, in which the flux qubits 
are placed inside two tunably coupled TLRs via a dc-SQUID. 

Our paper is organized as follows. In Sec.|ll] one-step gen- 
eration of the GHZ state in flux qubits coupled to a TLR sys- 
tem is presented using both analytical and numerical analysis. 
In SecHni one-step generation of the GHZ state in flux qubits 
coupled to two or more TLRs with tunable coupling is dis- 
cussed. Finally, discussion and conclusion are presented in 
SecHYl 



II. A GHZ STATE GENERATION FOR FLUX QUBITS 
INSIDE A CAVITY 



A. Theoretical Model 

The investigated system is schematically shown in Fig. [T] 
where superconducting flux qubits are strongly coupled to a 
one-dimensional transmission line resonator (TLR) with the 
geometric length Lq, distributed inductor L, and capacitance 
C. We consider the fundamental mode of the TLR, which can 
be modeled as a simple harmonic oscillator with the Hamilto- 





TLR 



(1) 



where we have set ?i = 1, the frequency of the fundamental 
mode is given by ujr/2n = 1/^/LC, and a are the cre- 
ation and annihilation operators of the fundamental mode of 
the resonator 

As shown in Fig.lifb), we assume that each flux qubit has 
a gradiometric configuration as studied in Ref. ^2^. The gra- 
diometric design is used to trap an odd number of fluxoids to 
bias the flux qubit at its optimal point, and also to greatly re- 
duce environmental noise. In contrast to the flux qubit with 
three junctions, one of which has a critical current a times 
smaller than that of the two identical junctions, here the small 
a-junction is replaced by a so-called a-loop, formed by a 
SQUID with two identical Josephson junctions. In this case, 
the ratio a and thus the coupling strength A between two 
circulating current states can be tuned via the magnetic flux 
/q$o through the SQUID loop, where $o is the magnetic flux 
quantum. The flux qubit Hamiltonian can be written as 



A(/a$o)o-x], 



(2) 



where the Pauli matrices read = |0)(0| — |1)(1| and = 
0) (1| + |1) (0|, in which |0) and |1) are the clockwise and coun- 
terclockwise persistent current states. fa) — 2/p/e$o is 
the biased magnetic energy where Ip is the persistent current 
in the qubit and A = /i — /2 is the magnetic frustration dif- 
ference in the two loop halves of the gradiometer. The energy 
gap A(/q,$o) can be controlled through the external magnetic 
flux /a$0- 

Let us now assume that the flux qubits are placed around the 
current anti-node of the TLR. The coupling between qubits 
and the TLR is approximately homogeneous because the di- 
mension of the qubits is on the order of several micrometers. 
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FIG. 1: Schematic description of our setup, (a) The flux qubits are 
placed around current anti-node of the TLR (half wavelength or full 
wavelength). The red dashed curve illustrates the amplitude of the 
magnetic field of the half-wavelength TLR. (b) Schematic descrip- 
tion of the tunable gap gradiometric flux qubits. The crosses denote 
Josephson junctions. The qubit is controlled by the bias lines and 
la- (c) Energy-level diagram of the system. A, lj^, i^r are energy 
gap of the qubits, driving microwave frequency, and fundamental fre- 
quency of the cavity mode, respectively. The strong classical driving 
field denoted by the large arrow resonantly interacts with the flux 
qubits. We assume that all qubits have the same Rabi frequency fJ/j. 
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which is much smaller than the wavelength of a few centime- 
ters for the fundamental electromagnetic modes in microwave 
frequency regime. We also assume that the distance between 
the two nearest flux qubits is sufficiently large 80 /im) 
to ensure that there are no direct interaction between dif- 
ferent qubits. The coupling strength between the k-th flux 
qubit and the TLR is given as = Mklpiro, where Mk is 
the mutual inductance between the qubit and the resonator, 
Ira = y/huJr/L is the zero-point current in the resonator 

The flux qubits are assumed to be near the optimal point 
(e « 0), and thus the total system Hamiltonian is as follows 
(refer to see Appendix A for detailed derivations): 



Hs — Hq + Hi + Hd, 



(3) 



with 



N 



TT t , - h 

fc=i 

N 

Hi ^^gk{a+ + a)a'^, 

k=l 
N 

Hd ^^flRCOs{oJdt) 



(4) 
(5) 
(6) 



k=l 



Here, Hq is the free Hamiltonian of the qubits and the cavity 
mode. Hi is the interaction Hamiltonian between the qubits 
and the cavity mode, and Hd describes the interaction between 
the qubits and classical driving microwave field. We assume 
that all qubits have the same frequency Ud of the driving mi- 
crowave field and the same Rabi frequency fi/j. It should be 
noted here that the homogenous coupling fin is induced by 
the classical field that is applied to the TLR (see Appendix 
A). 

In the basis of the eigenstates of the qubits and neglect- 
ing fast oscillating terms using the rotating-wave approxima- 
tion (RWA) if gk ^ ujr and flfj ^ ujd, the Hamiltonian in 
Eq. (|3]l takes the form 



^ A ^ 
Hi = w^a+a + ^ + X! 9k{a^<j'l + aa\) 

k=l 



N 



fc=l 

iujdi I ^k ^i^di\ 



(7) 



fc=i 



Because the energy gap A^, for each qubit can be tuned by 
the biased flux in the a loop, without loss of generality, we 
assume the energy gap of flux qubits are equal, i.e., Ak = A. 
Thus, in the rotating reference frame at the frequency Ud — A, 
the Hamiltonian in Eq. (|7]) is changed to 

~ ^ ^ O 

H2^Sa+a + Y, 9k{a+<J>l + aa^) + ^ ^« + a1), 

fe=l k=l 

(8) 

where we have defined the detuning 6 = uir — ujd > between 
the cavity field and the driving field, and the ladder operators 



(7+ = |~)(+| and (T_ = |+)(— I by using ground 1+) and first 
exited |— ) states of qubits. (If (5 < 0, it only means the transi- 
tion frequency of flux qubits is larger than the frequency of the 
fundamental cavity mode and also works well in our scheme.) 
Here, we assume that the circuit QED system works in small 
detuning regime, i.e., S > gk (not S ^ gk)- The third term of 
the Hamiltonian in Eq. ([8]l implies the free Hamiltonian of the 
dressed qubits by the driving field^". 

In the interaction picture with the unitary transformation 



U{t) = exp{-iHht) 
for the free Hamiltonian 



(9) 



N 



Hi, = 5a+a + Y,^{aX+at), 



k=l 



the interaction part of the Hamiltonian in Eq. ^ 



N 



Hint ^J^dki^^^'^- + 



fc=l 



becomes 

% - uHt)Hi,,tU{t) 



N 

Eft 

fc=i 



e 



+ olcT;, - 1(7 J e 



-iQ.Bt 



]} + H.c. 



(10) 



In the strong-driving regime, i.e., VLn ^ 5, gk, we can ne- 
glect fast oscillating terms and the Hamiltonian in Eq. ( fTOl l 
turns into 



E9k kr 



I t iSt~\ 

+ a'e I 



(11) 



The time evolution operator for the Hamiltonian in Eq. (fTTT i 
can be writte n-^ '-^^ as 



U{t) 



f[exp(- [ Blit)dBk{t)a',ai] (12) 

N N 



with 



(13) 



It is obvious that Bk (t) is a periodic function and vanishes 
at t = T„ — 2mT/d with the integer n. At these times, U{t) 
is independent of the variables of the cavity field and the flux 
qubits are decoupled from cavity field. We define 



Bl{t)dB,{t) 



9kg J 

4(5 



iS 



iSt 



(14) 
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At time r„, the time evolution operator takes the form 



N 



= exp —I 



N 

E 



2(52 



■ (15) 



For the convenience of the discussions, let us now as- 
sume that the qubits are equally coupled to the resonator, i.e., 
gk = (jj = g, and thus we assume = 7. It should be no- 
ticed that the inhomogenous coupling is not a significant prob- 
lem in our scheme. Actually, the homogenous coupling is not 
really required for this type of gate. The extra phase resulted 
from the inhomogenous coupling can be easily corrected by 
single qubit operations as mentioned in experiments^-^. If 
we adjust the parameters such that 7 = (1 + 2m)TT/8 with 
an arbitrary integer to, then the initially unentangled state 
\tp{0)) = '^k=i\+)k can be changed to the GHZ state with 
the unitary evolution [7(7), here |±) = (|0) ± \1))/V2. The 
parameter 7 is a geometric phase which will be addressed 
elsewhere. If the parameters are selected as n = 1, m = 0, 
5fc = 27r X 50 MHz, S = 2gk, then the GHZ state 



IGHZ) . -L 




N 



k + e 



■i7r(jV+l)/2 , 



(16) 



k=l 



is produced at the time T = 10 ns. We notice that there are 
two theoretical works on realization of controlled phase gate 
based on the 7 with superconducting qubits^^. 



B. Numerical Simulation 

To verify the validity of the approach proposed here, we 
now present numerical calculations. We simulate the dynam- 
ics of the system by using both its full Hamiltonian in Eq. Q 
without making any approximation, and the effective Hamil- 
tonians in Eqs. ( fTOb and ( fTTT i. For the convenience, we can 
rewrite the full Hamiltonian in Eq. Q, in the basis of the 
eigenstates of the qubits and the rotating reference frame of 
ujd, as 



N 



(17) 



fc=i 



N 



J2 gfc(a+4e*("'-+"'*)* + aa'le 



fc = l 



Let us assume that the cavity field is initially in the vacuum 
state |0)c, and the fc-th qubit is initially in the state \+)k = 
(|0) + 1 1))/ \/2, that is, the initial state of the whole system is 



N 



N 



1 



1^(0)) = 0$ \+)k\0)c = -= (|0)fe + \l)k) \0)c- (18) 
fc=i k=i V ^ 




FIG. 2: Time dependence of the fidelity of the generated entan- 
gled state in two flux qubits coupled to one TLR. The following 
parameters were used: ujr = 27rxlOGHz, A = 27r x 10.1 GHz, 
g ^ 2tv X 50 MHz, 5 = -2tv x 100 MHz. The target EPR state is 

(I -|- -|-) 4- i| ))/ \/2- (a) The blue solid, red solid, and red dashed 

curves are simulated using the effective Hamiltonians in Eq. dllt and 
in Eq. l ilOt and the full Hamiltonian in Eq. l il7t . respectively, with the 
optimized driving strength Q — 205. (b) Simulation using the full 
Hamiltonian in Eq. J17t for different driving strengths around the 
time t — 2n/5 = 10 ns. 



We can solve the Schrodinger equation to obtain the state 
\ip{t)) at any time by using the Hamiltonians in Eq. ( fTol i. 
Eq. ( fTTT l and Eq. ( [TtI i. Then we compare these states 
with the expected ideal GHZ state by using the fidelity 



F{t) = Tr[pGHZPg{t)l 



(19) 



where Pq{t) is the reduced density matrix of the N qubits and 
Pghz is the density matrix of the A^-qubit GHZ state. 

As an example, let us consider the interaction between two 
qubits and the cavity field with the initial state | -I- +)|0)c. 
In Fig. 12] the fidelities F{t) of different outcomes versus the 
evolution time t are plotted. In Fig.|2ja), we find that both the 
full Hamiltonian in Eq. ( fTTT l and the effective Hamiltonians in 
either Eq. ( fTOl i or Eq. (fTTT l well describe the dynamics of the 
two-qubits system. The expected entangled state (| + +) + i| — 
— ))/\/2 can be generated in t 10 ns, which is comparable 
with the single-qubit operation time. 
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To illustrate how the driving strength affects the fidelity of 
the expected GHZ state, we plot Fig.|2lb) by using the Hamil- 
tonian in Eq. ( fTTI i for different driving strengths. It can be 
clearly seen that the counter-rotating term of the driving field 
becomes important and reduces the fidelity when the driving 
strength is too strong. On the other hand, if the driving is too 
weak and the condition ftj^ ^ gk,S is not met any more, the 
fast oscillating terms in Eq. (fTOl i must be taken into account. 
Therefore, the driving field has to be optimized so that the 
maximum fidelity can be achieved. We provide a more de- 
tailed analysis of the effect of these fast oscillating terms on 
the fidelity of the generated state below. With the optimized 
driving strength, the fidelity is at best above 99.5%. It should 
be noted here that we did not make an adiabatic approxima- 
tion for the fast variable 57 in our simulation. Therefore, to 
prepare high fidelity GHZ states, the accuracy of control of 
the microwave pulse time should be determined by 2-K/Vtji 
in our scheme from analysis of the full system Hamiltonian. 
For example, to realize the preparation with a fidelity exceed- 
ing 90%, the precision of the microwave pulse time should 
be around 100 ps for the detuning 5 = —2tt x 100 MHz and 
VI n = 27r X 2 GHz in Fig.|2b). This is easily realized with a 
commercial pulse generator having a precision of 10 ps. 

In Ref. 24, to increase the effective coupling between 
qubits, they have to increase g (the coupling strength between 
qubits and resonator) and decrease the resonator frequency 
Ur- We argue that it is valid in experiments to make an adia- 
batic approximation for the resonator frequency uJr when g 
is approaching lo^ as in Fig. 2 of Ref. '24. Otherwise, to 
prepare high fidelity GHZ states in their scheme, the accu- 
racy of control of the dc current pulse time is determined by 
1-K jijJr, not by ^-KbJrl g^ ■ It means that the accuracy must be 
on the order of a few hundred ps for the selected parameters 
= 27r X 1 GHz and g = 27r x 144 MHz in Fig. 2 of Ref. M 
Such accuracy would present a major challenge for state-of- 
the-art dc current pulse technology. 



Vi{t) - 1-i / dt'H^jit') 
Jo 

-ff dt'dt"Hnjit')H„j{t") + --- .(21) 
Jo Jq 



We can treat U{t) as a constant during the integration because 
Hn{t) is oscillating much faster than the propagator. Then, 
we get 



Ui{t) 



9 



N 



2 N 



k=l 



9 



Y,{{l~cos{2nRt))U\t)alalU{t)} 



(22) 



N 



Near the time t = 2TTn/5, U{i) w exp(— i ^ Tc^cr^) and 
we obtain the fidelity 



F{t) ^ 1 - ^^^^/^^ (1 - cos(21}flt)), (23) 



where N is the number of qubits involved in the gate. It 
should be noted that the estimation of F{t) is obtained in the 
interaction picture. In the rotating frame and near the time 
t = 27rn/(5, the time evolution operator takes the form 



N 



N 



U{t) « exp(-z ^ nRtal/2) exp(-i ^ 7a^ai). (24) 

k k^j 



C. Nonideal Case 

In the derivation of effective Hamiltonian in Eq. ( fTOl i, we 
have neglected the following terms 

N 

fc=l 

- ((T^^+z(T^)e-'"^*] ae-'** + if.c.} (20) 

N . 

= H ^ [sin(17Ki)CT^ - cos(f7flt)a^] ae"*** + H.c. . 

k=l 

These terms could reduce the fidehty of the gate. However, 
we can use the spin-echo technique to eliminate the errors 
from the az terms. Below, we will only study the effect of 
the ay terms on the gate operation using the method described 
inRef.l26i 

In the interaction picture, the interaction Hamiltonian is 

Hnj{t) = U\t)Hn{t)U{t), and we have the propagator 
Ui {t) from the Dyson series. 



In the experiment, we can accurately control the duration 
of classical microwave field to fulfill the condition both 5t = 
2mT and ilj^t = An'-K in and n' are arbitrary integers) so that 

N 

the effect of the ay terms and exp(— i ^ Vljita^/2) vanishes. 

fe 

When N is not too large, e.g., N < flr/g, F{t) is mainly 
limited by the accuracy of control of the microwave pulse. On 
the other hand, if the accuracy of control of the microwave 
pulse is fixed, there is a polynomial decrease of the fidelity 
according to TV. 

Our scheme works in the strong-driving regime, i.e., fi/f ^ 
6, gk, so we can select the optimized driving strength Vlji — 
20(5 = 40(7 from the numerical simulation in Fig.|2|b) and use 
Eq. ( l23T l and Eq. (l24l i to roughly estimate the fidelity of the 
generated multi-particle GHZ state. In this case, F{t) ^ 1 
near the time t = 2tt/6 = 10 ns for N < 40. We can conser- 
vatively expect that our one-step proposal will work well for 
more than 20 qubits are placed at the current anti-node of the 
resonator by considering the inhomogenous coupling between 
qubits and resonator 
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III. SCALABLE CIRCUIT WITH TUNABLE TLRS AND 
TUNABLE NODES 

A. Scalable model 



coupling constant between two TRLs, where Iao, Ibo are the 
zero point current in cavity A and cavity B, respectively. It 
is possible to tune Mos by changing the penetrated flux 0e in 
the dc-SQUID loop, given by^ 



In our design, the half wavelength of the TLR for the funda- 
mental mode frequency cj^ = 27r x 10 GHz is around 6 mwM. 
If the distance between the two nearest flux qubits is around 80 
yum, we can only place about 6 flux qubits around the current 
anti-node position where the current variation is about 0.5%. 
Although the electric field is not completely zero around the 
anti-node, its effect on flux qubits is negligibly small. It is pos- 
sible to place more than 6 flux qubits in the half wavelength 
TLR because we can, in principle, design the coupling induc- 
tance so that the coupling constant is uniform in a wider 
range around the current anti-node. If there are extra phases 
arising from inhomogenous coupling between the qubits and 
TLR, we can correct them with additional single-qubit opera- 
tions. 

Here, we study the one-step generation of high-fidelity 
GHZ states for many qubits. We did not select a longer res- 
onator with a low frequency of the fundamental cavity mode 
to solve the problem of the scalabihty, for the following two 
reasons. First, the lower frequency of the fundamental cavity 
mode may lead to more operational errors resulting from the 
thermal excitations in the cavity. Second, the derived effective 
Hamiltonian in Eq. (fTTT i is valid in the strong-driving regime. 
If the driving strength in units of Rabi frequency approaches 
the energy gap of qubits, the counter-wave terms have to be 
taken into account, as shown in Fig.|2b). 

To solve the problem of the scalability, let us now use two 
coupled TLRs as an example to show how the multiparticle 
GHZ state can be generated via several TLRs in one step. As 
shown in Fig. [3] N qubits are placed into two cavities formed 
by two TRLs coupled by a symmetric dc-SQUID. We assume 
that the qubits labeled from 1 to [iV/2] interacted with TRL A, 
and other qubits labeled from [N/2] to N are coupled to TRL 
B, where means the maximum integer no more than 

N/2. Near the optimal point and in the basis of the flux qubit 
persistent current states, the total system Hamiltonian can be 
given as 



H' 



fe=i 



k=l 



N 



N 



i=[f]+i j=[f]+i 

N 

+J{(t)e)ia''b + ah'') + E ^RCOs{uJdt)a'^. (25) 
fc=i 

Here, the operators a(a^) and b{b^) are annihilation (creation) 
operators for the field in cavity A with the frequency cJq and 
cavity B with the frequency oj^, respectively, gk is the cou- 
pling constant between the fc-th qubit and the TRL A, and gj 
denotes the coupling of the j-th qubit and the TRL B. 

The parameter J{(j)e) = Mcs{(j>e)lAolBo is the inductive 



Meff(</.e) 



Lr 



,(26) 



2 + Pl cos{Itt — TT— ^) 

$0 



where / is an arbitrary integer, Ic is the critical current of 
the two identical Josephson junctions in the dc-SQUID, and 
screening parameter = 2TTLJc/^a < 1 to ensure that the 
coupler works in the nonhysteretic regime. We have ignored 
the direct inductive coupling between two resonators. More 
detailed discussions of the dc-SQUID coupler can be found in 
Refs.|33and38. 

We assume that all qubits are equally driven by the classi- 
cal field with the frequency w^, and the coupling between the 
qubits and the driving field is characterized by the constant 
il i^ . We rewrite the Hamiltonian in Eq. ( |25l ) in the basis of the 
qubit eigenstates and make the rotating wave approximation, 
thus the Hamiltonian in Eq. dZSl l becomes 



Hi 



L 2 J A I 2 J 

.,ata + ^_^^fc + ^g^(^fc„t+7j.c.) (27) 



fc=i 



k=l 



JV N 



+J{cbe){ah + ab^) + E ^(^e-'""* + H.c). 



fc=i 



Without loss of generality, in the following discussion we 
assume uja = i^b = which can be experimentally achieved. 
Because the TLR is a distributed element, we can easily de- 
sign its fundamental mode frequency in experiments. Now 
we introduce a canonical transformation to Eq. ( |27] ) via the 
operators 



P 



(28) 



Since the operators a and b describe different cavity fields in 
the different cavities, they satisfy the condition [a, b] = and 
then [P, Q] = 0. Substituting Eq. (|28]l into Eq. (|27]) and in the 
rotating reference frame of the resonant frequency uJd = A, 
we have 
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■HSHi — 

Qr cos ujdt Resonator A 



a 



out 



lfl+1 



M, 



CA 



CB 



N 



Coi 



ut 



Resonator B cos ujdt 



FIG. 3: Schematic description of flux qubits coupled to two tunably coupled TLRs (half-wavelength TLRs or full-wavelength TLRs). N qubits 
are coupled to the current anti-node of the two resonators. A dc-SQUID is placed at the current anti-node of the two resonators. The inductive 
coupling between two resonators could be tuned by changing the biased flux in the symmetric dc-SQUID loop. Two classical microwave 
fields resonantly interact with A'^ qubits by driving the TLRs. Lc, Mqa, Mqb are self-inductance of the dc-SQUID coupler loop, mutual 
inductance between the coupler and Resonator A, and mutual inductance between the coupler and Resonator B, respectively. 



[N/2] 



k=l 



j = [N/2] + l 

, lN/2] 

^{io - J)Q^Q + ^ E + 4Q) 

fc=l 



W N 
j = [W/2]+l 



k=\ 



(29) 



Using the same method as for the derivation of the effective 
Hamiltonian in Eq. (fTTT i for the flux qubits coupled to a TLR, 
based on Eq. ( |29l ). we can derive an effective Hamiltonian for 
the qubits coupled to two TRLs as 



[JV/2] 



N 

- E 

j=lN/2] + l 



V2 



-93 



k=\ 



(30) 

with 8' = uj — ujd- Here, we already neglect fast oscillating 
terms under the strong driving condition fi/j ^ gk, S' ,J. 



When S' > J, gk, and S' ~ ( £^ is an arbitrary odd inte- 
ger), at the time t = Tn ~ 2?™/ J, the flux qubits are decou- 
pled from resonators. In this case, we have 



^ 1 



[W/2] 

E 

N 

E 



(^' + J)(5'- J) 



gkgj 



{S' + J){5' - J) 



-5' 



[N/2] N 

E E 



fc=i 1+1 



{5' + J){5' -J) 



J 



(31) 



In the above equation, the first and second terms show that 
the flux qubits coupled to the same resonator are in permu- 
tation symmetry. The third term shows that the coupling be- 
tween any two qubits mediated by the two coupled TLRs is 
weaker than the first two terms, and there is also a sign dif- 
ference between the first two terms and the third term. If we 
adjust the detuning 5' and the coupling constant J{4>e) such 
that the conditions 

and 

are satisfied simultaneously, where I and m are arbitrary inte- 
gers, then the multiple particle GHZ state can be generated in 
one step. If we chose n = 1, m ~ I = Q, J = 2tt x 40 MHz, 
gk = V2J(k = 1 • • • iV), (5' = -3 J = -27r x 120 MHz, the 
GHZ state can generated in 25 ns. We would like to emphasize 
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Resonator I 




Resonator M 



FIG. 4: Schematic description of a circle array of TLRs coupled by 
a dc-SQUID. A/ TLRs are placed as the current anti-nodes around 
the coupler. Flux qubits are coupled to the current anti-nodes of the 
TLRs. 



that J = 27r X 40 MHz is a reasonable value in experiments. 
If we assume Iayq = Ibio = 50 nA, J — 2ttx40 MHz means 
that the effective mutual inductance Mcs between two TLRs 
mediated by dc-SQUID is around 5.32 pH. It is experimen- 
tally realizable with the selected parameters as Ic = 1.5 fi A, 
McA = McB = 60 pH and Lc = 200 pH for the dc-SQUID 
coupler It is possible to continuously tune the coupling from 
anti-ferromagnetic to ferromagnetic like a rf-SQUID coupler 
in experiments^^. We notice that there is a theoretical work 
on a flux qubit mediating the coupling between two TLRs and 
working as a quantum switch^. In that scheme, it is also pos- 
sible to realize tunable coupling between two TLRs mediated 
by the flux qubit. 

Now, we consider how to generate A^-qubit (N > 20) GHZ 
states in one step. As shown in Fig. |4l a dc-SQUID coupler 
is coupled to M TLRs which form a circle array. Flux qubits 
are placed at current anti-nodes of the TLRs. For the conve- 
nience of discussions, we assume that the coupling constants 
between resonators and the center coupler are homogenous. 
The resonators are in permutation symmetry. If we select pa- 
rameters to simultaneously meet conditions such as those in 
Eq. d32] i and Eq. ( [33T l, we can, in principle, extend our scheme 
to the case of many qubits. 



B. Numerical Simulation 

For flux qubits coupled to a system of two coupled TLRs, 
the numerical simulation procedure is similar to that in pre- 
vious case. We can assume that all qubits are initially in the 
state <E)^^i\+)k and that the two cavity fields are initially in 
the vacuum state |00)c, i.e., that the initial state of the whole 
system is |-0(O)) = «)f=il+)fc|00)c. The full Hamiltonian 




t(ns) 

FIG. 5: Time dependence of the fidelity of the generated entan- 
gled state for two flux qubits which are coupled by two coupled 
TLRs. Here, the two qubits are separately coupled to the two TRLs. 
We set the following parameters for our numerical calculations: 
oj^ = 27r X 10 GHz, A = 2% x 10.12 GHz, J = 27r x 40 MHz, 

Qk = \/2Jand5' = -3J. The target state is + })/V2. 

(a) The blue solid and red dashed curves are simulated using the ef- 
fective Hamiltonian Eq. l l30b and the full Hamiltonian in Eq. i34\ . 
respectively, with optimized driving strength Qr = 42 J. (b) Sim- 
ulation using the full Hamiltonian in Eq. J34l ( with different driving 
strengths around the time t = 2n/ J = 25 ns. 

used in the simulation is 

~ ^ o 

HUs^m = + + (34) 

k=l 
, lN/2] 

fe=i 

] = lN/2]+l 
, lN/2] 

3 = {Nl2\+\ 
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Comparing the simulation results shown in Fig. |5]with the 
full Hamiltonian and the effective Hamiltonian, we find that 
the effective Hamiltonian in Eq. dSOl l can also describe the dy- 
namics of two qubits system well. With the optimized driving 
strength, the maximum fidelity can be above 99.8%. The two 

qubits entangled state (| + +) + j| ))/V2 can be produced 

mt — 2TT/J — 25 ns. Our proposal provides an obvious ad- 
vantage in that the multiparticle GHZ states are generated in 
one step in ^ 25 ns with high fidelity. In contrast, in the ex- 
periments described in Refs. 14 and 15, three-qubit GHZ state 
was generated step by step and the total generation time was 
^ 80 ns with a fidelity of ^ 90%. According to those experi- 
ments, the total generation time will linearly increase with the 
numbers of qubits. 

IV. DISCUSSION AND SUMMARY 

Let us now discuss the experimental feasibility. For a TLR 
with the fundamental mode frequency = Stt x 10 GHz 
and the quality factor Q ^ W^, the decay rate is about ^ 
1 MHz. The decoherence time achieved in the tunable gap 
gradiometric qubit is longer than 1 /is in experiments"*'. Thus 
the GHZ state can be generated in about 10 ns, which is much 
shorter than the decoherence time of qubits and the decay time 
in the TLR. 

To demonstrate the GHZ state, we could use the well devel- 
oped quantum state tomography technique in superconducting 
qubits**^ to reconstruct the density matrix of the final state'^''*-'. 
After the GHZ state is generated, we can tune the transition 
frequencies of the qubits such that the interactions between 
all qubits and the TLRs are switched off. For example, the 
change of detuning from 2n x 100 MHz up to 27r x 2 GHz 
reduces the coupling strength by a factor of 20. It is suffi- 
cient for experimental demonstration of the GHZ state using 
quantum state tomography. 

We derive the effective Hamiltonian under the strong- 
driving condition. The strong driving on the flux qubit as high 
as 27r X 2 GHz has been realized in our group. A detailed 
theoretical analysis of leakage to higher energy levels under 
strongly resonant microwave driving is given in Ref. 44. From 
these calculation, we can conclude that the leakage to higher 
energy levels is negligibly small because of the strong anhar- 
monicity for flux qubits at the optimal point (e « 0). 

In summary, we have proposed a scheme for generation of 
the multiqubit GHZ entangled state in tunable gap flux qubits 
coupled to a TLR. We also extend this scheme to the case of 
tunable flux qubits coupled to two or more coupled TLRs. 
The operation time is comparable to that of the single-qubit 
operation. In principle, our proposal can be generalized to the 
case in which the tunable qubits are coupled to a circle array 
formed by AI coupled TLRs via a dc-SQUID coupler, that is, 
our proposal is scalable to some extent. We expect that our 
scheme has an useful contribution to the generation of cluster 
states for one-way quantum computation. 



Appendix A: Homogenous coupling between qubits and 
classical driving field through TLR 

Let us now show how to obtain the Hamiltonian in Eq. (O 
when the classical driving field is applied to the TRL. In our 
assumption, the driving microwave with large amplitude is 
substantially detuned from the frequency of resonator In this 
situation, quantum fluctuation in the drive is very small com- 
pared to the drive amplitude, and the drive field could be con- 
sidered as a classical field"*^. The classical microwave field 
driving on the resonator can be described by 

Hd = u{t)a+e-"^-'' + V* (i)ae*"''* (Al) 

where v{t) is the amplitude and ijjd is the frequency of exter- 
nal driving. The Hamiltonian of the whole system is H — 
Ho + Hi + Hd, where Hq in Eq. (gl and Hj in Eq. Q. We 
can then displace the field operators using time dependent dis- 
placement operator 

£i(/3) = exp(/3a+ - l3*a), (A2) 

where /3 is an arbitrary complex number and the field a goes 
to a + /? under this unitary transformation. 

In the case in which the driving amplitude e is independent 
of time, the displaced Hamiltonian in the energy eigenstate 
basis of qubits reads 

H = D+{l3)HD{l3)~tD+{/3)D{fi) (A3) 

^ N N 

fc=i fc=i 

^ N N 

= a;^a+a + -^Afe(T^ + ^5fc(a++a)cr^ + 

k=l k=l 

2gkV cos{uJdt) f, 
-5 

fe=i 

where /3(t) = ve^^'^'^^/S, 5 = ojr - oJd- 

The above Hamiltonian is the Hamiltonian in Eq. ([3]) 
in energy eigenstates basis of the qubits with ~ 
2gkV cos{ujdt) / 5. If the qubits are equally coupled to the res- 
onator, i.e., gk = g, we obtain homogenous coupling be- 
tween qubits and classical driving field through the TLR. 
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